Homodyne tomography characterization and nonlocahty of a dual-mode optical qubit 
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A single photon, delocalized over two optical modes, is characterized by means of quantum ho- 
modyne tomography. The reconstructed four-dimensional density matrix extends over the entire 
Hilbert space and thus reveals, for the first time, complete information about the dual-rail optical 
quantum bit as a state of the electromagnetic field. The experimental data violate the Bell inequality 
albeit with a loophole similar to the detection loophole in photon counting experiments. 
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According to Feynman, single-particle interference is 
"a phenomenon which is impossible, absolutely impossi- 
ble, to explain in any classical way^and which has it in 
the heart of quantum mechanics" The explanation 
offered by quantum mechanics is that a particle incident 
onto a beam splitter is not, as classically expected, ran- 
domly reflected or transmitted, but forms a delocalized 
coherent superposition of these two possibilities. The 
quantum state in the two beam splitter output modes A 
and B can be expressed as 



\^)=t\Ia,^b)-p\^aAb) 
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where and are, respectively, the beam splitter trans- 
mission and reflectivity, and the kets in the right-hand 
side are written in the particle-number basis. 

In this article, we investigate the delocalized state 
formed by a photon. We employ homodyne tomog- 
raphy, a state characterization technique based on 
phase-sensitive measurements of the electromagnetic field 
quadratures. Although there is only one light particle, we 
find it to affect the field in both modes simultaneously, 
giving rise to nonclassically correlated, phase-dependent 
quadrature statistics. This is a direct consequence of the 
coherent, entangled nature of the state and is remark- 
able because the optical phase in both the single-photon 
and vacuum states considered individually is completely 
uncertain. 

We show that the delocalization of the photon in the 
state can be formulated as a noncontextual violation 
of local realism, namely as a violation of Bell's inequality 
Nonlocahty of the single photon has been discussed 
earlier in a number of theoretical publications, and vari- 
ous experiments were proposed for its demonstration 
none of which have so far been realized. The present 
work achieves this goal by using a very different measure- 
ment method and by introducing an assumption which 
resembles the fair sampling assumption in traditional ex- 
periments on nonlocahty 4]. To our knowledge, this is 
the first experimental verification of Bell's theorem per- 
formed in a continuous- variable setting. 

Apart from its fundamental implications, our experi- 
ment finds its use in the linear optical implementation of 
quantum optical information processing [3 , in which the 



state (0 plays the role of a quantum bit. To date, char- 
acterization of optical qubits has been based on studying 
relative photon number statistics in each mode and in 
their various linear superpositions as well as (in the case 
of multiple qubits) photon number correlations between 
modes. Employing this approach. White et al. have im- 
plemented tomography of entangled two-qubit systems 
6]. Most recently, this technique has been extended 
to characterization of quantum dynamical processes 0. 
Similar methods are used in celebrated quantum cryp- 
tography and quantum teleportation protocols. 

A major drawback of the photon counting approach to 
quantum state characterization is the a priori assumption 
that the modes involved are either in one of the states 
|1a,0b), |0a,1b) or in their linear combination. As a 
result, being characterized are not the quantum states of 
the carrier modes, but their projections onto a subspace 
spanned by these basis vectors. It is neglected that in 
actual experimental situations, most of the time the state 
of the modes does not belong to this subspace but is the 
double vacuum |0a,0b). These overwhelming events, as 
well as all other events not falling into the qubit frame, 
were simply eliminated from the analysis. 

The homodyne tomography method employed in this 
paper is free from this disadvantage as it permits com- 
plete characterization of a quantum optical ensemble. So 
far it has been employed to study single- 8] and multi- 
mode squeezed states and, more recently, the single- 
photon and related states Here we apply this tech- 
nique to the delocalized photon state (QJ and evaluate 
its density matrix and its Wigner function. A detailed 
theoretical analysis of different aspects of an experiment 
such as ours was made by Jacobs and Knight 11] as well 
as Grice and Walmsley 12]. 

In our experiment, the initial single-photon state was 
prepared by means of a conditional measurement on a 
biphoton produced via parametric down- conversion |l^ . 
We used frequency-doubled 2-ps pulses from a mode- 
locked Ti: Sapphire laser running at A = 790 nm which 
underwent down-conversion in a BBO crystal, in a type-I 
frequency-degenerate configuration |l(]| . 

Field quadratures of the delocalized photon ensem- 
ble were measured by means of two homodyne detec- 



tors (associated with fictitious observers Alice and Bob) 
placed into each beam splitter output channel (Fig. 1(a)). 
The local oscillator fields were provided by the master 
Ti: Sapphire laser and contained several million photons 
per pulse. Their optical mode had to be matched to those 
of the state under investigation . A specially designed 
balanced detector employed two Hamamatsu S3883 pho- 
todiodes of a 94-% quantum efficiency whose photocur- 
rents were directly subtracted |0. The amplification 
circuit provided suppression of the electronic noise in a 
large frequency range from DC to 2-3 MHz, which en- 
abled time-resolved quadrature measurements for each 
pulse. The time-integrated homodyne detector output 
directly corresponded to a quadrature value, with a scal- 
ing reference obtained from the vacuum state shot noise. 

With every incoming photon, both detectors made a 
measurement of the field quadrature Xa and Xb with 
the local oscillators' phases set to Oa and 9b^ respectively. 
The quadrature statistics collected at various phases were 
used to reconstruct the density matrix of the generated 
ensemble. 

Generally, homodyne tomography reconstruction of a 
two-mode state requires fine control of both local oscilla- 
tor phases involved. However, the particular state stud- 
ied in the present experiment was generated by splitting 
a single photon which has no optical phase. The sum of 
Alice's and Bob's phases Oa-\-Ob is therefore meaningless 
and does not affect the homodyne statistics. Therefore, 
we let this phase vary randomly and controlled only the 
relative phase SO = Oa — Ob of the two modes. During an 
experimental run, we varied this parameter slowly over 
a 27r range, and acquired a few hundred thousand pairs 
(Xa^Xb) from the homodyne detectors. We have exe- 
cuted two data acquisition runs using two different beam 
splitters with transmissions equal to 0.5 and 0.08. 

Fig. l(c,d) shows histograms of the dual- mode quadra- 
ture measurements. The two-dimensional distribution 
Pt§0{Xa, Xb) indicates the probability of detecting a 
particular pair {Xa^Xb) of quadratures at a given local 
oscillator phase setting. These densities are the marginal 
distributions of the four-dimensional Wigner function of 
the two-mode ensemble being measured. They can be 
used for its direct reconstruction via the inverse Radon 
transformation. 

In this work, however, he have used a more precise 
maximum-likelihood reconstruction technique ITS ll^ . 
We applied the iterative algorithm described in to the 
two-mode density matrix of the reconstructed ensemble 
in the Fock representation: 
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with the photon number restricted to five per mode. 
The positive operator- valued measure (POVM) operator, 
used to describe the homodyne detector, was modified 
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FIG. 1: (a) Scheme of the experimental setup. LOs: local 
oscillators, (b-d) Histograms of the experimental quadra- 
ture statistics pi^q(Xa, Xb) for the "zero reflectivity" (b), 
symmetric (c), and highly reflective (92%, d) beam splitters. 
Phase-dependent quadrature correlations are a consequence 
of the entangled nature of the state |^). Also shown are in- 
dividual histograms of the data measured by Alice and Bob, 
which are phase-independent. 



with respect to that of an ideal homodyne detector in two 
ways. First, we corrected for the non-perfect detection 
efficiency associated with linear optical losses and non- 
unity quantum efficiency of the homodyne detector pho- 
todiodes. These losses amounted to 1 — 77det = 0-14 and 
are described mathematically in the form of a Bernoulli 
transformation affecting every photon entering the detec- 
tor 17, 18, 19]. 

Secondly, due to a random variation of the phase 
Oa^Ob^ all the POVM elements, for which m-\-n ^ k-\-l, 
vanish, and so do the respective elements of the recon- 
structed density matrix — akin to the off-diagonal terms 
of any phase-randomized single-mode ensemble. 

The reconstructed density matrix (Fig. 2) features 
a strong contribution of the double-vacuum term 
|0,0)(0,0|. This is a consequence of imperfect prepara- 
tion of the initial single photon: instead of the state |1), 
a statistical mixture = r] |1)(1| + (1 — ^) |0)(0| is avail- 
able at the beam splitter input |l(i[l4]. The vacuum frac- 
tion is directly transferred to the dual-mode ensemble: 
p = + (1 - 7?) |0,0)(0,0|. The reconstructed en- 

semble is in excellent agreement with the above equation 
and corresponds to the preparation efficiency r] = 0.64. 

Some insight into the quadrature dynamics of the non- 
local single-photon state can be gained by analyzing its 
Wigner function (WF). Theoretically, the transformation 
of a quantum ensemble by a beam splitter can be ex- 
pressed as a rotation of its WF in the four-dimensional 
phase space 0: 

W{Xa, Pa.Xb, Pb) = WUtXa + pXb, (3) 
tPa + pPb.tXb - pXa, tPb - pPa)- 



3 




FIG. 2: Density matrix (absolute values) of the measured 
ensemble for the symmetric (a) and highly reflective (b) beam 
splitters in the photon number representation (|^. 




FIG. 3: Cross-sections of the two- mode Wigner function 
W{Xa, Xb, Paj Pb) of the reconstructed ensemble for the 
symmetric (left column) and highly reflective (right column) 
beam splitters, (a, b): cutting plane (Xa = Pa = 0); the WF 
is axially symmetric, (c, d): cross-sections through the plane 
(Pa — Pb — 0) possesses mirror symmetry with symmetry 
planes oriented at an angle tan~^(y9/r) with respect to the 
coordinate axes. Insets: cross-sections through Xb = 0. 



In our experiment, the input WF is an uncorrelated prod- 
uct of the Wigner functions of the single-photon and vac- 
uum states in the beam splitter input modes: 

W^n (Xi , Pi , X2 , P2 ) = 1 ) (Xi , Pi ) X 1^1 0) (X2 , P2 ) . (4) 

A rotation in the phase space entangles the input mode, 
generating the ensemble |^). The effect of this rotation 
can be seen in Fig. 3 which shows the cross-sections of 
the four-dimensional WF of the reconstructed ensemble 
through two different planes. As expected for the sin- 
gle photon, the WF is negative in the phase space ori- 
gin, albeit the negativity is reduced due to a non-perfect 
preparation efficiency. 

A similar argument can be used to understand the 
shape of the measured quadrature probability densities. 
Consider an experiment in which the beam splitter has 
zero reflectivity, so that Alice measures the field quadra- 
tures of the single-photon state, and Bob, at the same 
time, those of the vacuum state. The joint quadra- 
ture statistics are, naturally, uncorrelated [Fig. 1(b)]. 
However, if the two modes are subjected to a beam 
splitter transformation prior to their measurement, their 
field amplitudes will mix, the quadrature histogram will 
rotate by an angle tan~^(r//)) and become correlated 
[Fig. l(c,d)]: 

Xa = rXi+pX2 (5) 
Xb = -pXi+rX2. 

The above relation is valid only if Alice and Bob mea- 
sure the same quadrature of the incoming fields, i.e. for 
SO = OT 60 = TT. On the other hand, the setting 



50 = 7t/2 corresponds to a measurement of orthogo- 
nal quadratures, e.g. X by Alice and P by Bob. It 
can be interpreted as a simultaneous, unprecise measure- 
ment of the canonically conjugated position and momen- 
tum quadratures of the input single-photon state with 
an eight-port homodyne detector [l^,!!^. Furthermore, 
the distribution pr^/2(^A, Xb) acquired with a symmet- 
ric beam splitter [Fi g;. 1 (c), middle] is equal to the Q- 
function of this state [21j . 

As discussed above, these phase-dependent features of 
the quadrature statistics are incompatible with a naive 
picture according to which the photon is randomly local- 
ized in one of the beam splitter output channels. The 
nonlocal character of our measurements can be formal- 
ized as a violation of the Bell inequality 0. In order to 
apply the Bell theorem to the continuous field quadra- 
tures measured in this experiment, the latter had to be 
converted to a dichotomic format. We proceeded as fol- 
lows. Suppose the output of each homodyne detector is 
processed by a discriminator producing the value S = 1 
if the observed quadrature exceeds a certain threshold 
T, = — 1 if it is below — T, and generating no out- 
put otherwise. Restricting to the events in which both 
discriminators have generated a value, we studied corre- 
lations Eab =< SaSb > between the values acquired by 
Alice and Bob. 

The correlation Eab as a function of Oa — Ob exhibits 
quantum interference [Fig. 4] which is similar to that ob- 
served in Bell's Gedankenexperiment for the spin cor- 
relation as a function of the angle between the axes of 
the Stern-Gerlach apparata. Increasing the threshold T 
leads to a higher visibility of the interference pattern. For 
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FIG. 4: Bell inequality violation by a threshold-selected sub- 
set of the experimental data [symmetric beam splitter], (a) 
Interference pattern for the threshold of T — 0.85 shows an 
amplitude of 0.818 di 0.018; (b) theoretical and experimen- 
tal dependence of the amplitude on the threshold value. A 
higher threshold increases the amplitude, but also reduces the 
number of selected data. 



T > 0.54, the amplitude exceeds 1/V2 and the Bell in- 
equality is violated. Interpreting our data in this manner, 
we observed violations by up to six standard deviations. 

One might object that such an interpretation involves 
rejection of some of the the acquired data. This objection 
is valid; however, similar selection is also a part of the 
traditional way of verifying Bell's theorem by studying 
the counting statistics of a polarization-entangled pho- 
ton pair Q. Indeed, every single-photon counter has a 
discriminator as a part of its electronic circuit, and the 
signals not exceeding the discrimination threshold are 
disregarded. Viewed outside of a context of any partic- 
ular physical theory, both tests possess an equal degree 
of validity: they both involve a fair sampling assumption 
and thus suffer from the detection loophole [2^ . 

Is our entire data set noncontextually nonlocal? Al- 
though a number of schemes of reformulating Bell's the- 
orem for continuous variables (2^ have been proposed, 
we have not yet found one that could be applied to our 
original data. On the other hand, we could not come up 
with a local hidden variable model that would replicate 
the quadrature statistics observed in this experiment 24] . 

In summary, we have characterized a two-mode op- 
tical qubit using homodyne tomography, reconstructing 
its density matrix and the Wigner function. For the first 
time, complete information about this quantum ensem- 
ble is revealed, including those terms that do not belong 
to the qubit subspace of the Hilbert space. The experi- 
mental data demonstrate a nonlocal character of the de- 
localized single photon but the question of noncontextual 
nature of this nonlocality remains open. 
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